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Abstract 



We revisit the application of different separability criteria by recourse to 
an exhaustive Monte Carlo exploration involving the pertinent state-space of 
pure and mixed states. The corresponding chain of implications of different 
criteria is in such a way numerically elucidated. We also quantify, for a 
bipartite system of arbitrary dimension, the proportion of states p that can 
be distilled according to a definite criterion. Our work can be regarded as a 
complement to the recent review paper by B. Terhal [Theor. Comp. Sci. 287 
(2002) 313]. Some questions posed there receive an answer here. 
Pacs: 03.67.-a; 89.70.+c; 03.65.Bz 
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I. INTRODUCTION 



The development of criteria for entanglement and separability is one aspect of the current 
research efforts in quantum information theory that is receiving, and certainly deserves, 
considerable attention [1]. Indeed, much progress has recently been made in consolidating 
such a cornerstone of the theory of quantum entanglement [1]. The relevant state-space 
here is of a high dimensionality, already 15 dimensions in the simplest instance of two-qubit 
systems. The systematic exploration of these spaces can provide us with valuable insight 
into some of the theoretical questions extant. 

As a matter of fact, important steps have been recently made towards a systematic 
exploration of the space of arbitrary (pure or mixed) states of composite quantum systems 
[5-7] in order to determine the typical features exhibited by these states with regards to the 
phenomenon of quantum entanglement [5-11]. Entanglement is one of the most fundamental 
and non-classical features exhibited by quantum systems [12], that lies at the basis of some 
of the most important processes studied by quantum information theory [12-16] such as 
quantum cryptographic key distribution [17], quantum teleportation [18], superdense coding 
[19], and quantum computation [20,21]. 

It is well known [1] that, for a composite quantum system, a state described by the density 
matrix p is called "entangled" if it can not be represented as a mixture of factorizable pure 
states. Otherwise, the state is called separable. The above definition is physically meaningful 
because entangled states (unlike separable states) cannot be prepared locally by acting on 
each subsystem individually [22]. 

The separability question has quite interesting echoes in information theory and its asso- 
ciate information measures or entropies. When one deals with a classical composite system 
described by a suitable probability distribution defined over the concomitant phase space, 
the entropy of any of its subsystems is always equal or smaller than the entropy character- 
izing the whole system. This is also the case for separable states of a composite quantum 
system [23,24]. In contrast, a subsystem of a quantum system described by an entangled 
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state may have an entropy greater than the entropy of the whole system. Indeed, the von 
Neumann entropy of either of the subsystems of a bipartite quantum system described (as 
a whole) by a pure state provides a natural measure of the amount of entanglement of 
such state. Thus, a pure state (which has vanishing entropy) is entangled if and only if its 
subsystems have an entropy larger than the one associated with the system as a whole. 

Regrettably enough, the situation is more complex when the composite system is de- 
scribed by a mixed state: there are entangled mixed states such that the entropy of the 
complete system is smaller than the entropy of one of its subsystems. Alas, entangled mixed 
states such that the entropy of the system as a whole is larger than the entropy of either of 
its subsystems do exist as well. Consequently, the classical inequalities relating the entropy 
of the whole system with the entropy of its subsystems provide only necessary, but not 
sufficient, conditions for quantum separability. There are several entropic (or information) 
measures that can be used in order to implement these criteria for separability. Consider- 
able attention has been paid, in this regard, to the g-entropies [1,24-31], which incorporate 
both Renyi's [32] and Tsallis' [33-35] families of information measures as special instances 
(both admitting, in turn. Shannon's measure as the particular case associated with the limit 
g — > 1). The reader is referred to Appendix A for a brief review on g-entropies. 

The early motivation for the studies reported in [24-31] was the development of practical 
separability criteria for density matrices. The discovery by Peres of the partial transpose 
criteria, which for two-qubits and qubit-qutrit systems turned out to be both necessary and 
sufficient, rendered that original motivation somewhat outmoded. In point of fact, it is not 
possible to find a necessary and sufficient criterion for separability based solely upon the 
eigenvalue spectra of the three density matrices pab,Pa — TtbIpab], and ps — TtaIpab] 
associated with a composite system A(B B [23]. 

Interesting concepts that revolve around the separability issue have been developed over 
the years. A beautiful account is given in Terhal in [1]. Among them we find criteria like 
the so-called Majorization, Reduction and Positive Partial Transpose (PPT) together with 
the concept of distillability [1] . Quantum entanglement is a fundamental aspect of quantum 
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physics that deserves to be investigated in full detail from all possible points of view. The 
chain of implications, and the related inclusion relation, among the different separability 
criteria is certainly a vantage point worth of detailed scrutiny. It is our purpose here to 
revisit, with such a goal in mind, the separability question by means of an exhaustive Monte 
Carlo exploration involving the whole space of pure and mixed states. Such an effort should 
shed some light on the inclusion issues that interest us here. Concrete numerical evidence 
will thus be provided on the relations among the separability criteria. We will then be able 
to quantify, for a bipartite system of arbitrary dimension, the proportion of states p that can 
be distilled according to a definite criterion. This numerical exploration could be viewed as 
a complement on the review paper by Terhal [1] , because some questions posed by her will 
receive an answer in this work. 

The paper is organized as follows. We sketch in Section II the different separability 
criteria to be investigated and discuss some mathematical and numerical techniques used in 
our survey in Section III. Our results are reported in Section IV, and some conclusions are 
drawn in Section V. For the sake of completeness, we include an Appendix on g-entropies. 

II. BRIEF SKETCH ON SEPARABILITY CRITERIA 

Prom a historic viewpoint, the first separability criterion is that of Bell (sec [1] and 
references therein). For every pure entangled state there is a Bell inequality that is violated. 
It is not known, however, whether in the case of many entangled mixed states, violations 
exist. There does exist a witness for every entangled state though [4]. It was shown by 
Horodecki et al. that a density matrix p = pab is entangled if and only if there exists an 
entanglement witness (a hermitian super-operator W — W^) such that 

TrWp<0, while 

TrW p > 0, for all separable states. (1) 
A special, but quite important LOCC operational separability criterion, necessary but 
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not sufficient, is provided by the positive partial transpose (PPT) one. Let T stand for 
matrix transposition. The PPT requires that 

[i®f](p)>0. (2) 

Another operational criterion is called the reduction criterion, that is satisfied, for a given 
state p = pab, when both [1] 

i ® - P > 

PA^i-p>0. (3) 

Intuitively, the distillable entanglement is the maximum asymptotic yield of singleton 
states that can be obtained, via LOCC, from a given mixed state. Horodecki et al. [36] 
demonstrated that any entangled mixed state of two qubits can be distilled to obtain the 
singleton. This is not true in general. There are entangled mixed states of two qutrits, 
for instance, that cannot be distilled, so that they are useless for quantum communication. 
In our scenario an important fact is that all states that violate the reduction criterion are 
distiUable [37]. 

Entanglement witnesses completely characterize the set of separable states. Alas, they 
are not usually associated to a simple computational treatment, except in the PPT instance. 
Thus, in order to decide whether a given state p is entangled one needs additional criteria, 
functional separability ones [1]. One of them associates PPT to the rank of a matrix. 
Consider two subsystems A, B whose description is made, respectively, in the Hilbert spaces 
Tin and Tim- Focus attention now in the density matrix p = pab for the associated composite 
system. If 

1. p has PPT, and 

2. its rank is such that TZ < max(n, m), 

then, as was proved in [38], p is separable. The above referred to entropic criteria are also 
functional separability ones. Still another one is majorization. 
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Let {Aj} be the set of eigenvalues of the matrix and {7^} be the set of eigenvalues of 

— * 

the matrix ^2- We assert that the ordered set of eigenvalues A of ^1 majorizes the ordered 
set of eigenvalues 7 of ^2 (and writes A :^ 7) when J2^=i Xi > J2i=i li for all k. It has been 
shown [23] that, for all separable states pab = P, 

>- Ap, and 

Kb ^ (4) 

There is an intimate relation between this majorization criterion and entropic inequalities, 
as discussed in [1,24]. 

III. SEPARABILITY PROBABILITIES: EXPLORING THE WHOLE STATE 

SPACE 

We promised in the Introduction to perform a systematic numerical survey of the prop- 
erties of arbitrary (pure and mixed) states of a given quantum system by recourse to an 
exhaustive exploration of the concomitant state-space S. To such an end it is necessary to 
introduce an appropriate measure // on this space. Such a measure is needed to compute 
volumes within S, as well as to determine what is to be understood by a uniform distribu- 
tion of states on S. The natm^al measure that we are going to adopt here is taken from the 
work of Zyczkowski et al. [5,6]. An arbitrary (pure or mixed) state /? of a quantum system 
described by an A^-dimensional Hilbert space can always be expressed as the product of 
three matrices, 

p - UD[{\}]Ul (5) 

Here U is an N x N unitary matrix and D[{Aj}] is an N x N diagonal matrix whose diagonal 
elements are {Ai, . . . , Ajv}, with < Aj < 1, and J2i \ — 1. The group of unitary matrices 
U{N) is endowed with a unique, uniform measure: the Haar measure v [39]. On the other 
hand, the TV-simplex A, consisting of all the real /^-uples {Ai, . . . , Aat} appearing in (5), is a 
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subset of a {N — l)-dimensional hyperplane of TZ'^ . Consequently, the standard normalized 
Lebesgue measure C^-i on TZ^"^ provides a natural measure for A. The aforementioned 
measures on U {N) and A lead then to a natural measure /j, on the set S of all the states of 
our quantum system [5,6,39], namely, 

H^ux Cn-1- (6) 

All our present considerations are based on the assumption that the uniform distribution 
of states of a quantum system is the one determined by the measure (6). Thus, in our 
numerical computations we are going to randomly generate states according to the measure 
(6). 

IV. SURVEY'S RESULTS 

A. The overall scenario 

An overall picture of the situation we encounter is sketched in Fig. 1, that is to be 
compared to Fig. 3 of [1] . Notice that our numerical exploration allows us to dispense with 
TerhaVs interrogation signs. This constitutes part of the original content of the present 
Communication. 

The set of all mixed states presents an onion-hke shape, as conjectured by Terhal [1]. 
Which among these states are separable? As reviewed above, several criteria are available. 
We start with the g-entropic one (see the Appendix and [40]). By using a definite value of q, 
namely g = oo, and the sign of the associated, conditional g-entropy, we are able to define a 
closed sub-region, whose states are supposedly separable. This region has a definite border, 
that separates it from the sub-region of states entangled according to this criterion. What 
we see now is that, if we use now ot/ier separability-criteria, the associated sub- regions shrink 
in a manner prescribed by the particular criterion one employs. The shrinking process ends 
when one reaches the sub-region defined by the Positive Partial Transpose (PPT) criterion. 
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which is a necessary and sufficient separabihty condition for 2 x 2 and 2x3 systems, being 
only necessary for higher dimensions. 

Summing up, the volume of states which are separable according to different criteria 
diminish as we use stronger and stronger criteria. There is a first shrinking stage associated 
to entropic criteria, from its Von Neumann (g = 1) size, as q grows, to the limit case q ^ oo 
[40] . A second stage involves majorization, reduction, and, finally positive partial transpose 
(PPT) [1]. 

B. PPT and Reduction 

We report now on our state-space exploration with regards to the probability of finding 
a state with positive partial transpose. The results are depicted in Fig. 2. The solid fine 
corresponds to states with dimension N = 2 x N2, while the dashed fine corresponds to 
N = 3 X N2 states. Note how similar are the pertinent values in both cases. The tiny 
difference between them can be inspected in the inset (a semi- logarithmic plot). To a good 
approximation, our PPT probabilities decrease exponentially. 

Fig. 3 deals instead with the probability of finding a state which obeys the strictures of 
the reduction criterion, ior N — 2xN2 (sohd fine) and N — 3xN2 (dashed line) . As a matter 
of fact, PPT and reduction coincide ior N — 2x N2. It is known that if a state satisfies PPT, 
it automatically verifies the reduction criterion [1]. Here we have demonstrated that, at least 
in the N — 2 x A^2-instance, the converse is also true. However, in the TV = 3 x A^2-case, it 
is much more likely to encounter a state that verifies reduction than one that verifies PPT. 

C. Entropic criteria and Majorization 

We begin with a brief recapitulation of former g-entropic results. The situation encoun- 
tered in [41] was that the "best" result within the framework of the "classical g-entropic 
inequalities" as a separability criterion was reached using the limit case q ^ 00, but con- 
siderably less attention was paid to other values of q. This was remedied in [40], where 
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the question of g-entropic inequalities for finite g- values was extensively discussed. It was 
there re-confirmed that the above mentioned limit case does indeed the better job as far as 
separability questions are concerned [40]. For such a reason, this limit g- value is the only 
one to be employed below. See the Appendix for more details on g-entropies. 

In Fig. 4 we depict the probability of finding a state which, for g — > oo, has its two 
relative g-entropies positive (dashed curves). In view of the intimate relation of entropic 
inequalities with majorization [1,24], we also analyze in Fig. 4 the probability that a state 
is completely majorized by both of their subsystems (sohd line). It is shown in [24] that, if 
Pab satisfies the reduction criterion, its two associated relative g- entropies are non-negative 
as well. 

In the same work the authors assert that majorization is not implied by the relative 
entropy criteria. Our results confirm this assessment. In Fig. 4, the lower curves correspond 
to states p with N — 2 x N2, while the upper curves have N — 3 x N2. Majorization 
results and g-entropic do coincide for two-qubits systems (A^i — N2 — 2). More generally, 
majorization probabilities are a lower bound to probabilities for relative g-entropic positivity, 
an interesting new result, as far as we know. Notice also that the two approaches yield quite 
similar results in the N — 3 x N2 case. 

D. Comparing more than two criteria together 

We compare now the reduction criterion to the PPT one. The former is implied by the 
latter but is nonetheless a significant condition since its violation implies the possibility of 
recovering entanglement by distillation, which is as yet unclear for states that violate PPT 
[24] . Fig. 5 a) depicts the probability that state p with N — 3 x N2 either: 

1. has a positive partial transpose and does not violate the reduction criterion, or 

2. has a non positive partial transpose and violates reduction. 
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Remember that in the case N — 2 x N2, the two criteria always coincide [1]. For 3 x N2 the 
agreement between the two criteria becomes better and better as A^2 augments. 

Of more interest is to compare the relations among PPT, majorization, and the entropic 
criteria (Fig. 5b), since it is not yet known how the majorization criterion is related to other 
separabihty criteria hke PPT, undistillability, and reduction [24]. In this vein. Fig. 5 b) 
plots the "coincidence-probability" between, respectively, 

1. PPT and majorization (solid line), and 

2. PPT and the q'-entropic criterion (dashed line). 

The curves on the top correspond to N — 2 x N2, while those at the bottom to A?^ = 3 x A^2- 
In this last case the two curves agree with each other quite well. 

The conclusion here is that, as N2 augments, the probability of coincidence among the 
three criteria, and in particular between majorization and PPT (our main concern), rapidly 
diminishes at first, and stabilizes itself afterwards. For two qubits the three criteria do agree 
with each other to a large extent. 

Fig. 6 a) depicts the probability that, for a given state p, 

1. reduction and majorization (solid line) and 

2. reduction and the g-entropic criterion (dashed line) 

yield the same conclusion as regards separability. Without PPT in the game, and opposite 
to what we encountered in Fig. 5, we find better coincidence ior N — 3 x N2 systems (top) 
than for N = 2 x N2 (bottom). The deterioration of the degree of agreement as N2 grows is 
similar to that of Fig. 5, though. 

Fig. 6 b) represents the probability that a state, for q — > 00, either: 

1. has both positive relative g-entropies and satisfies the majorization criterion, or 

2. has a negative relative g-entropy and is majorized by both of their subsystems. 
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The solid line corresponds to the case N — 2 x N2, while the dashed lines corresponds to 
the N — 3 X N2 instance. These results together with those of Figs. 4-5 could be read as 
implying that majorization and the g-entropic criteria provide almost the same answer for 
dimensions greater or equal than N — 3 x N2. 

Finally, in Fig. 7 we look for the probability Pagree that all criteria considered in the 
present work do lead to the same conclusion on the separability issue. Pagree is plotted as a 
function of the total dimension N — Ni x N2, with Ni — 2 (solid hne) and Ni — 3 (dashed 
line). The agreement is quite good for two qubits, deteriorates first as A^2 grows, and rapidly 
stabilizes itself around a value of 0.26 for Ni = 2 and of 0.1 for A^i = 3. 

E. Distilling 

Let us at now consider the results plotted in Fig. 8. We ask first for the relative number 
of states that violate the reduction criterion and are thus distillable [36] (solid line), and 
appreciate the fact that, as N grows, so does the probabihty of finding distillable states. On 
the other hand, the probability of encountering states that violate the majorization criterion, 
represented by dashed lines, is much lower than that associated to distillation. 

For both criteria, the upper solid line corresponds to the case N — 2 x N2, and the lower 
one to N — 3 X N2. The dashed curve with crosses represents the case A'" = 2x N2, while the 
one with squares indicates the N = 3x N2 instance. The dependence with N2 of the dashed 
curves (majorization violation) is not so strong as that of the solid ones (distillability) . Our 
results are lower bounds to the total volume of states that can be destilled. 

V. CONCLUSIONS 

We have performed a systematic numerical survey of the space of pure and mixed states 
of bipartite systems of dimension 2x N2 and 3 x 7V2 in order to investigate the relationships 
ensuing among different separability criteria. Our main results are 
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• Regarding the line of separability implication, see our graph in Fig. 1 and compare 
with the similar one of Terhal's (her Fig. 3). Her interrogation signs receive an answer 
in our Fig.l 

• It is known that if a state satisfies PPT, it automatically verifies the reduction criterion 
[1]. In the present work we show that in the N — 2x 7V2-instance, the converse is also 
true. In the A^" = 3 x A^a-case, it is much more hkely to encounter a state that verifies 
reduction than one that verifies PPT. 

• We have numerically verified the assertion made in [24] that majorization is not implied 
by the relative entropic criteria. Majorization results and g-entropic criteria coincide 
for two-qubits systems. In general, majorization probabilities constitutes lower bound 
for relative g-entropic positivity. 

• Regarding the relation between majorization and PPT, the agreement between the 
criteria deteriorates as N2 grows. 

• For dimensions > 3 x iV2, as illustrated by Figs. 4-5, majorization and the g-entropic 
criteria provide almost the same answers. 

The present authors believe that the results of this numerical exploration shed some light 
on the intricacies of the separability issue. 
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FIGURE CAPTIONS 



Fig. 1 - Schematics of the inclusion relations among separbility criteria as given by the 
volume occupied by states p for a given dimension A^" which fulfill them. 

Fig. 2 - Probability of finding a state with positive partial transpose. The solid line corre- 
sponds to states with dimension N — 2xN2, while the dashed fine corresponds to A'" = 3 x A^2 
states. The difference between these curves can be appreciated in the inset (semi-logarithmic 
plot). Our probabilities decrease, to a good approximation, in exponential fashion. 

Fig. 3 - Probability of finding a state fulfilling the reduction criterion for N = 2 x N2 (solid 
line) and N — 3 x N2 (dashed fine). The two probabihties coincide ior N — 2 x N2. 

Fig. 4 - Probability of finding a state whose two relative g-entropies are positive for g — > 00 
(dashed curves). The probability that a state be completely majorized by both of their 
subsystems is represented by the solid line. Bottom: curves correspond to states p with 
N^2x N2. Top: N^Sx N2. 

Fig. 5 a) Probabihty that the state p with N — 3 x N2 either has i) a positive partial 
transpose and does not violate the reduction criterion, or ii) has a non positive partial 
transpose and violates reduction. In the case N — 2 x N2 the outcome is always unity. Fig. 
5 b) Probability that 

1. PPT and majorization (sohd line) and, 

2. PPT and the g-entropic criterion (dashed line) 

lead to the same conclusion regarding separabihty. Top: N — 2 x N2. Bottom: N — 3 x N2. 

Fig. 6 a) Probability that reduction and majorization (solid line) and reduction and the 
g-entropic criterion (dashed line) yield the same conclusion reagrding separability. Top: 
N = 3 X N2. Bottom: N = 2 x N2 (lower curves). Fig. 6 b) Probability that a state, for 
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g — > oo, either i) has both positive relative g-entropies and fulfills majorization, or ii) has 
a negative relative g-entropy and is majorized by both of their subsystems. The solid line 
corresponds to the case N = 2 x N2, while the dashed line corresponds to N = 3 x N2. 

Fig. 7 - Total probability that all criteria considered in the present work lead to the same 
conclusion regarding separability. Probabilities are plotted as a function of the total dimen- 
sion N ^ NiX N2, with iVi = 2 (solid fine) and A^i = 3 (dashed line). 

Fig. 8 - Solid line: probability that a state violates the reduction criterion. Dashed line: the 
same for violation of the majorization criterion. Top: N — 2 x N2. Bottom: N — 3 x N2. 
The dashed curve with crosses represents the case N — 2 x N2, while the one with squares 
indicates the N = 3 x N2 instance. 

APPENDIX A: Q-INFORMATION MEASURES AND THE ISSUE OF 

QUANTUM SEPARABILITY 

There are several useful entropic (or information) measures for the investigation of a 
quite important subject: the violation of classical entropic inequalities by quantum entangled 
states. The von Neumann measure 



is important because of its relationship with the thermodynamic entropy. The g-entropy, 
which is a function of the quantity 



provides one with a whole family of entropic measures. In the limit g — > 1 these measures 
incorporate (Al) as a particular instance. Most of the applications of g-entropies to physics 
involve either the Renyi entropy [32], 



Si = -Tr(plnp), 



(Al) 



= Tr (p^) 



(A2) 




1 



In (Ujg) , 



(A3) 



1-g 
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or the Tsallis entropy [33-35] 

5^ = ^(1-^.)- (A4) 

We reiterate that the von Neumann measure (Al) constitutes a particular instance of both 
Renyi's and Tsallis' entropies, obtained in the limit g — > 1. The most distinctive single 
property of Tsallis' entropy is its nonextensivity. The Tsallis entropy of a composite system 
A® B whose state is described by a factorizable density matrix, p^B — Pa® Pb, is given by 
TsaUis' g-additivity law, 

Sf\pAB) = S^^\pa) + Sf\pB) + (1 - q)S^rMSf\pB). (A5) 

In contrast, Renyi's entropy is extensive. That is, if pab — Pa® Pb, 

5f (PA^) = 5f (pa) + Sf\pB). (A6) 

TsaUis' and Renyi's measures are related through 

5f ) = ), (A7) 

where the function F is given by 

F{x) = {e('-«)^ - l} . (A8) 

An immediate consequence of equations (A7-A8) is that, for all non vanishing values of q, 
Tsallis' measure S*^"^-* is a monotonic increasing function of Renyi's measure S^^\ 

Considerably attention has been recently paid to a relative entropic measure based upon 
TsaUis' functional defined as 



1 + (1 -gjS"^ (Pb) 



Here pab designs an arbitrary quantum state of the composite system A®B, not necessarily 
factorizable nor separable, and pb = TrAipAs)- The relative g-entropy Sj^^\B\A) is defined 
in a similar way as (A9), replacing pB by pA — Ttb^Pab)- The relative g-entropy (A9) 
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has been recently studied in connection with the separability of density matrices describing 
composite quantum systems [29,30]. For separable states, we have [24] 

5f (A|5)>0, 

S^^\B\A)>0. (AlO) 

On the contrary, there are entangled states that have negative relative g-entropies. That is, 
for some entangled states one (or both) of the inequahties (AlO) are not verified. 
Notice that the denominator in (A9), 

1 + (1 - g)5f ) ^ wg > 0. (All) 

is always positive. Consequently, as far as the sign of the relative entropy is concerned, 
the denominator in (A9) can be ignored. Besides, since Tsallis' entropy is a monotonous 
increasing function of Renyi's (see Equations (A7-A8)), it is plain that (A9) has always the 
same sign as 

5f {A\ B) = Sf> (pab) - 5f ) (pb). (A12) 



19 



figl 



states p 



^-Entropic 




4 8 12 16 20 




i 



4 8 12 16 20 

N 




4 6 8 10 12 14 16 18 20 

N 



